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I Introduction 



A Historical Background 

The concept of differential invariants was first introduced by B. Riemann in 1858 in his 
classical works^'^ concerning the Euler equations for an ideal fluid flow in two independent 
variables 

p'(p) 

ut + uux H px = 0, pt + up^ + pw^ = 0, p > 0. (1) 

P 

Here, u is the local velocity of the fluid, the pressure p is assumed to be a different iable 
function of the density p and p' = dp/ dp. Riemann investigated the formulation and 
mathematical correctness of problems involving the propagation and superposition of 
waves described by equation ([T]). Next, he constructed rank-2 solutions corresponding to 

1 /2 

the "superposition" of two waves propagating with local velocity u = ± {p'{p)) when 
the directions of propagation are parallel and opposite. The main element of Riemann's 
method is the introduction of new dependent variables (called Riemann invariants) which 
have the property of preserving their values along appropriate characteristic curves of the 
original system. This approach enables us to reduce the number of dependent variables 
for the wave superposition problem. Consequently it simplifies the task of solving the 
original system. In the case of Euler's equations ([1]) this method allows us to construct a 
scattering double wave solution 

ii = A;i/2(ri-r2)+Mo, p = Ae''+'\ p = kAe'-'^'' + po, uo,Po,k,AeR, (2) 

which reduces ([1]) to the following invariant hydrodynamic system 

rl + {k^/\r^-r^ + l)+Uo)rl = 0, r^^ + {k^^^r^ - - 1) + Uq) rl = 0. (3) 

Riemann studied the asymptotic behavior of initial localized disturbances (i.e. initial 
data with compact support) corresponding to the above-mentioned waves. Consequently, 
he demonstrated that even for sufficiently small initial data, after some finite time T, 
these waves could be separated again in such a way that waves of the same type as 
those assumed in the initial data could be observed. Riemann noticed that the solution 
for systems of hydrodynamic type ([3]), even with arbitrarily smooth initial data, usually 
could not be continued indefinitely in time. After a certain finite time T the solutions blew 
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up (more precisely, the first derivatives of the considered solution become unbounded after 
some finite time T > 0). So, for times t > T, smooth solutions to the Cauchy problem 
do not exist. This phenomenon is known as the gradient catastrophe.^'"^ Furthermore, 
Riemann was interested in extending the solution in some more generalized sense beyond 
the time T of the blow up. On the basis of the conservation laws for the mass, energy 
and momentum, Riemann^ and later Hugoniot^ introduced the concept of weak solutions 
(non-continuous) in the form of shock waves. Based on this concept he proved some laws 
connecting the wave front velocity and parameters of the fluid state before and behind 
that discontinuity. 

The problem of propagation and superposition of Riemann waves has been investigated 
by many authors since then (see e.g. Ref.^"^^). A number of attempts to generalize the 
Riemann invariant method and its applications can be found in the recent literature on 
the subject (see e.g. Ref.^^"^^ and references therein). 

B Supersymmetric models 

Over the last three and a half decades, there has been a great deal of interest in super- 
symmetry. Supersymmetric theories involving both bosonic (even) and fcrmionic (odd) 
degrees of freedom have been used extensively to describe various types of physical phe- 
nomena. The construction of supersymmetric extensions of existing field theories, an 
approach first employed in the context of particle physics, has been successfully used to 
supersymmetrize physical theories involving both classical and quantum fields.^^'^^'^^ In 
particular, these techniques have been applied extensively to the area of fiuid mechanics, 
starting with the study of simple Euler-type equations.^^'^°'^^ A supersymmetric exten- 
sion of the Korteweg-dc Vrics equation was formulated through the use of a superspace 
and superfield mechanism,^^'^^ and a version of the symmetry reduction method adapted 
to Grassmann-valued partial differential equations was used to determine its invariant 
solutions. More recently, supersymmetric generalizations of the Chaplygin gas were for- 
mulated in (1 + 1) and (2 + 1) dimensions by supplementing the Lagrangian with additional 
fermionic fields. ^^'^^'^^ This latter approach was also used to build an = 1 supersym- 
metric extension of polytropic gas dynamics,^'' a covariant and supersymmetric theory of 
relativistic hydrodynamics in four-dimensional Minkowski space,^®'^^ and a Kaluza-Klein 
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model of a relativistic fluid. ^° 

In the last couple of years, there has been much interest in systems of partial differential 
equations describing a steady, irrotational and compressible fluid flow in a plane. '^^ Such 
a system can be written in the form 

Uy-Vx = 0, 

(4) 

{pu)x + {pv)y = 0, 

where [u, v) are the Cartesian components of the fluid velocity expressed in terms of a 
velocity potential ip (where u = ipx, v = ipy), and the density p is deflned as a function of 
u and V. Two particular cases stand out in physical importance. In the case where the 
density is given by 

P = e-", (5) 

the fluid equations (jll) correspond to those of a Gaussian, irrotational, compressible fluid 
flow 

(l - {(fxf) (fixx - 2!fx<fy<fxy + (l - ify?) fyy = 0. (6) 

Likewise, in the case where the density is 

p=ii+u'+vr'^', (7) 

the fluid equations (jlj) reduce to the minimal surfaces equation in (2 + l)-dimensional 
Minkowski space^^ 

^ + f I =0 



:i + + {f.ry'y , V(i + + w)"v, 

which can be linked through the Wick rotation y = it, to the scalar Born-Infeld equation 

(l + (iPxY) ^tt - 2ipxipt(Pxt - (l - i^tY) ^xx = 0. (9) 

In both cases (Gaussian fluid flow and scalar Born-Infeld equations ([2])), N = 1 
supersymmetric generalizations were formulated^^'^^ through the use of a superspace and 
superfleld formalism similar to that used for the Korteweg-de Vries equation. 

It should be noted that the Born-Infeld equation (Q is compatible with the following 
hydrodynamic type system expressed in terms of Riemann invariants'^ 

Rt + SRx = 0, St + RSx = 0. (10) 
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the two being linked through the transformation 

(1 + (^^2 _ (^^)2)l/2 ^^^^ 



R 



^ _ {I + {^.? - {^t?y/' ^.^t ^ ' 

Numerous multiple wave solutions of system ffTOl) can be determined from the invariant 
solutions of the Born-lnfeld equation through the transformation (ITT!) . However, since 
the explicit expressions of such solutions are complicated, they are not presented here. 

The system fllUp admits infinitely many conservation laws of the form 

^ + JC^) = 0, A;GZ+, (12) 
where the densities p^^-* and fluxes J*-'^-' are given by 

k /k-1 \ 

= R's''-\ f ''^ = -RS R'S''-'-' . (13) 
(=1 \j=i / 

Equations ([T2|) and (fT3|) imply the existence of closed 1-forms dx^''^ on f2 C X 

= -jWrft + pWrfx. (14) 

From the closure of the 1-forms dx^''^ (i-e. d{dx^''^) = 0), it follows that the integrals 

{-j'-'^Ut + p^'^^dx) =x'-''\x,t), (15) 



locally depends only on the end points of the curve 7 (i.e. it is locally independent of the 
trajectory in the plane R^). Each integral defines a mapping 

: (x,t)^x^'\x,t)eR. (16) 

We treat the real- valued functions x^'^'' i^, the coordinates of a two-dimensional surface 
JF immersed in M^. The map x = {x^^\ ■ ■ ■ iX^^^) is called the Weierstrass formula for 
immersion. 

It was shown^^ that all rank-2 solutions (i.e. the general integral) of the system (11 01) 
can be obtained by solving the following implicit system with respect to the variables R 
and S 

x = RF[{R)-F,{R) + SF^{S)-F2{S), t = F[{R) + F^{S), (17) 

where Fi and F2 are arbitrary differentiable single-variable functions of R and S, re- 
spectively. Here, the symbol prime means differentiation with respect to the argument 



{R or S). Hence, from the functionally independent functions R and S satisfying equa- 
tions f|T7|) we can construct formally, using the Weierstrass formula for immersion f|T5|) . a 
two-dimensional surface in M'^. 

C Objectives 

The primary objective of this paper is to construct a supersymmetric extension of the 
evolutionary system ( fTOl) involving Riemann invariants. For this purpose, we make use of 
a superspace and superfield formalism in which the space of independent and dependent 
variables is enlarged so that it also includes fermonic Grassmann variables. We seek to 
identify the symmetries, subalgebra classification and invariant solutions of the classical 
and supersymmetric models, and to compare the classical and supersymmetric cases. This 
work is a follow-up of research performed in Ref.^^ 

This paper is organized as follows. In Section II, we identify the symmetry algebra, 
subalgebra classification and solutions of the classical form of equation (flOll . In Section III, 
we present the proposed supersymmetric extension for the Riemann invariant system ffTOl) . 
In Section IV, we describe the superalgebra of Lie symmetries and perform a systematic 
classification of its one-dimensional subalgebras. In Section V, we use the symmetry 
reduction method to obtain various classes of invariant solutions of our supersymmetric 
model. Finally, Section VI contains a summary of our results and possible future outlook. 
Here, the symmetries and solutions of the classical Riemann system are compared to those 
of its supersymmetric extension. 

II Symmetry properties and solutions of the classical system in 
Riemann invariants 

In section, we describe the Lie algebra of symmetries of the classical version of the hy- 
drodynamical model ffTOl) and perform a systematic classification of its one-dimensional 
subalgebras. For each representative subalgebra, the symmetry reduction method is used 
to transform the system ffTOj) to a reduced system of ordinary differential equations. This 
allows us to generate solutions of (fTOl) which are invariant under the action of the partic- 
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ular subalgebra under consideration. Further details on the symmetry reduction method 
can be found in Ref.^^'^^ 

If we restrict our considerations to continuous Lie groups and use an infinitesimal 
approach, we find that the classical Lie point symmetry algebra L of the system fllOl) is 
spanned by the following six vetor fields 

ri = 9,., To = dt, W = td, + dR + ds, Mi = xd, + tdt, 

(18) 

M2 = xd^ - tdt + 2RdR + 2Sds, J = xdt - R^Sr - S^ds, 
where = d/dx, etc. The physical interpretation of this Lie algebra as it applies to 
the space of coordinates (x, t, R, S) is as follows. The vector fields Ti and Tq generate 
translations in the x and t directions respectively, the elements Mi and M2 are associated 
with two different types of dilations and W represents a modified Galilean boost. The 
generator J corresponds to a type of inverse boost where the roles of the space and time 
coordinates are exchanged. The commutation relations for the generators (1181) of the Lie 
algebra L are given in Table I. 

In order to perform a classification of the subalgebras of L, we make use of the following 
decomposition 

L = {{Ml} © {M2, W,J}^3 {Ti, To}, (19) 
where the symbol © represents a direct sum and 5 a semi-direct sum. 

We now proceed to classify the subalgebras of the Lie algebra L. That is, we construct 
a list of representative subalgebras of L such that each subalgebra of L is conjugate to 
one and only one element of the list under the equivalence conjugation relation 

X e^Xe-^ = X + [F, X] + i [Y, [Y, X]] + i [Y, [Y, [Y, X]]] + . . . (20) 

The reduction of the classical hydrodynamical system (|TOl) to systems of ordinary differ- 
ential equations requires consideration of the subalgebras of L. In order to perform the 
classification, we made use of the techniques for direct and semi-direct sums of algebras 
as described in Ref.^^ and references therein. For the sake of brevity, only the results are 
listed. The splitting one-dimensional subalgebras of L are 

Li = {Ti}, L2 = {Mi}, L3 = {M2}, L, = {W}, L, = {W-J}, 

L6 = {M2 + A;Mi, A; ^0,1,-1}, L7 = {Mi + M2}, Lg = {M2 - MJ, (21) 

Lg = {W + eMue = ±1}, L^q = {W - J + kMi,k ^ 0}, 

7 



and the non-splitting one- dimensional subalgebras of £ are 



Lu = {W + eTi}, Lu = {Mi + M2 + eTo}, = {M2 - Mi + eT^} . (22) 



We now proceed to use the symmetry reduction method in order to obtain invariant 
solutions of the classical hydrodynamical system (fTOj) . First, we find for each of the 
subalgebras listed in this section the associated three invariants along with the change of 
variable that must be substituted into system (ITUI) in order to obtain the set of reduced 
ordinary differential equations. In each case, the invariant involving only the independent 
variables (called the symmetry variable) is labelled by the symbol a. The invariants and 
change of variable are listed in Table II, while the system of reduced ordinary differential 
equations are listed in Table III. Where it is possible, we also give a solution of the reduced 
system in closed form. 

For subalgebra Li, we obtain the trivial solution where R and S are constants. This 
solution is also present in the case of subalgebra L2. The other solution which satisfies 
the reduced equations for L2 is the one where both R and S are equal to the simple ratio 
x/t, which represents a centered wave. 

The solution of the reduced equations in the case of subalgebra L4 is the following 
time-modified centered wave 

R(x,t) = - + Cit+—, S{x,t) = --Cit + —, (23) 
t ij t t 

where Ci and C2 are arbitrary constants. For subalgebra Lu, a modified version of this 
expression is obtained 

R{x, t) = -^ + Ciit + e) + -^, Six, t) = -^- Ci{t + e) + ^. (24) 

Exponential solutions are present in the case of subalgebras L-j and Lg, given by 

_ Cia:e^^(*+^^) _ CiX 

R{x,t) - ^Ciit+C2) _ 1' S{x,t) - ^(jiit+Ci) _ 1' 



and 



R{x, t) = ^ (e'^^^^^^^) - 1) , Six, t) = ^(l- e-^^(^+^^)) , (26) 



respectively. 
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For subalgebra L3, we obtain the following solution expressed in terms of radicals 

R{x, t) =1 + (^-^C2 - iCi|) (^-\c,xt + IC2 + y^^^^T^^^^^t+^Tie^^ , 

S(x, = - 4 ( ((4 - CiC2)xt + Ci)^CfxH^ - 2CiC2xt + C| + 16xt 



+ ((Ci^Ca - ACi)xH^ + (I2C2 - 2CiCi)xt + C|) ) x 



(^{-Cfxt + C1C2 - A)^ClxH^ - 2CiC2xt + C| + 16xt 
+ {CfxH^ + (12Ci - 2CfC2)xt + {CiCi - 4C2)) 



(27) 



For the remaining subalgebras, we did not obtain explicit solutions in closed form. We 
must therefore describe our results separately for each case. 

First, for subalgebra L5, the solution to the reduced equations is given by 

(sinF \ 

where F{a) is solved implicitly by the equation 

(1 - {smFf)^K^ - (sinF)2 + (sinF)^ = i(Co(7 + + 1) (29) 

Next, in the case of subalgebra Lq, the solution is given by 

= - ~. (30) 



where F{a) is the solution of the equation 

In the case of subalgebra Lg, the solution to the reduced equations is 

G{a)^a + e-^, (32) 

where F{a) satisfies the equation 

{sF - ea - l)F„a + {F - a){F„f = (33) 



For subalgebra Lio, the solution to the reduced equations is given by 

= arctan$((j), G{o-) = arctanr((T), (34) 

where 

= (35) 
and $(o") is solved implicitly by the differential equation 

-(T^(A; + $)(l + $2)$^, + (T^(2A;$ + 2$2-fc2-l)($,)2 + (T(A:-$)(l + $2)<l><,-(l + $2)2 = q 



(36) 



In the case of subalgebra L12, the solution is given by 



G{a) = (37) 
where F{a) is solved implicitly by the equation 

aF{eF - 2)F,, + 2a{F^y + 2F{eF - 1)F, = (38) 



Finally, for subalgebra L13, the solution is 

G{a) = 1^ ( ^ 
where F{a) satisfies the equation 



aF^ 



(39) 



aF{l + 2eF)F^^ - 2(t(1 + eF){F„f + 2F(1 + eF)F„ = (40) 

Unfortunately, the ordinary differential equations (l33l) , (!36|) , (!38|) and (HOj) do not have 
elementary solutions since they do not have the Painleve property. Thus, it is not possible 
to obtain solutions in closed form for the corresponding subalgebras. 



Ill = 1 supersymmetric version of the system in Riemann 
invariants 

Our purpose in this section is to supersymmetrize the hydrodynamic type system (fTOl) . 
where R{x,t) and S{x,t) are the standard (bosonic) Riemann invariant fields. In order 
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to accomplish this, we enlarge the space of independent variables {x, t} to a superspace 
{x,t,6}, where 6 is an independent odd (fermionic) Grassmann variable. In addition, we 
introduce the two fermionic superfields 

^x,t,e) = ^{x,t) + eR{x,t), '^{x,t,e) = ^/j{x,t) + es{x,t), (4i) 

where ^(x, t) and ip{x,t) are two new fermionic- valued fields. 

We want to build our supersymmetric extension in such a way that it is invariant 
under the supersymmetry transformation 

X ^ x + r]0, 9^9-7], (42) 

which is generated by the vector field 

Q = 9d,-de. (43) 

It is also convenient to define the covariant derivative 

D = 9d, + de, (44) 

which possesses the useful property that it anticommutes with the supersymmetry gener- 
ator Q. Thus, any system of superequations written in terms of the covariant derivative 
operator D and the superfields $ and \E' will be manifestly invariant under the supersym- 
metry transformation (H2ll . 

The most general form of the extension can be written under the form 

<^t + aD^D^$ + (1 - a)^£'^$ = 0, 

(45) 

+ bD^D^^ + (1 - b)^D'^^ = 0, 

where a and b are arbitrary real-valued constant parameters. When decomposed into 
coefficients of the various powers of 9, the system (145|) is equivalent to the following 
system of four differential equations expressed in terms of the four fields R{x,t), S{x,t), 
^{x, t) and i/j{x, t) 

Rt + SR:^ + aV^^^^ + (a - l)^'^^^ = 0, 
St + RS, + fe^x^x + l)e^xx = 0, 

(46) 

it + aSi^, + (1 - a)^R, = 0, 
ijt + bR^j, + (1 - b)iS, = 0. 
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In this paper, we consider the case where a = 1 and 6=1. Here, our supersymmetric 
hydrodynamical systems (H5|) and (l46l) become 

+ = 0, 

(47) 

+ D^D^-^ = 0, 

and 

Rt + SR, + ^^^x = 0, 

St + RS, + e.v^. = 0, 

(48) 

6 + = 0, 
ijt + Rip, = 0, 

respectively. In the limiting case where the fermionic fields ^ and ip approach zero, we 
recover the classical version of the system in Riemann invariants ffTOl) . 

IV Symmetries and subalgebra classifications 

In order to construct and investigate group-invariant solutions of the hydrodynamical 
supersymmetric system fHHj) by means of the symmetry reduction method, we need to 
find its Lie algebra C of infinitesimal symmetries and then classify all subalgebras of C of 
dimension 1 into conjugacy classes. 

A Symmetries 

The Lie superalgebra C of infinitesimal symmetries of the supersymmetric system fl48p is 
spanned by the following eight vector fields 

Po = dt, Pi = d,, Yi = d^, Y2 = d^, B = td, + dR + ds 

Di = 2xd, + 3tdt - Rdn - Sds, D2 = xd., + tdt + ^d^, (49) 

D3 = xd^ + tdt + ijjd^, 

where Pq, Pi, Yi and I2 represent translations in x, t, ^ and ip respectively, i? is a 
modified Galilean boost and Di, D2 and D3 correspond to three independent dilations 
in the independent and dependent variables. The commutation (and anticommutation in 
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the case of two fermionics) relations of the generators described in f H9|) are summarized 
in Table IV. 

The superalgebra L can be decomposed into the following composite semidirect sum 
C = l^[{Di,D2,Ds} B {fi}} 5 {Po,Pi}| 5 {Yi,Y2}. (50) 
It should be noted that this superalgebra is solvable. 

B Classification 

The one-dimensional subalgebras of the Lie superalgebra C can be classified in a manner 
similar to that used in Section II for the Lie algebra L of the classical hydrodynamical 
system. Again, we only list the results. The splitting one-dimensional subalgebras of C 
are 

£i = {Pi}, C2 = {Po}, C^ = {B}, C^ = {Di}, C, = {D2}, 
£e = {Ds}, Cj = {D2 + aDi,a^0}, Cs = {D^ + aDi,a 0}, 
£9 = {D3 + aD2, a ^ 0}, £10 = {D-s + aD2 + bDi,a,b^ 0}, 

(51) 

£11 = {D2 + eB,e = ±1}, £12 = {D3 + eB,e = ±1}, 
£13 = {D3 + aD2 + eB,a^O,€ = ±1}, £14 = {B + ePq, e = ±1}, 
£15 = {Fi}, £i6 = {1^2}, £17 = {Yi + eY2, e = ±1}. 
The non-splitting one-dimensional subalgebras of £ are 

£i8 = {Pi + mY, + rj2Y2}, £19 = {Po + ViYi + rj^2}. 

£20 = {B + rj^i + r^Fa}, £21 = {^1 + rj^i + ^^2], £22 = {D2 + rfe^a}, 

£23 = {^3 + mYi], £24 = {D2 + aDi + 772^2, a ^ 0}, (52) 

As = {^3 + aD^ + r^Fi, a ^ 0}, £26 = {D2 + sB + rf^2, e = ±l}, 

£27 = [D-i + eB + rnY,, e = ±1}, £28 = {B + ePo + rj^i + r]^2, e = ±1}. 

This classification allows us to obtain, systematically, all symmetry reductions of the 
hydrodynamical supersymmetric system fHHl) and to provide (where it is possible) solu- 
tions of the reduced equations. 
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V Invariant Solutions 



We now use the symmetry reduction method in order to obtain invariant solutions of 
the hydrodynamical supersymmetric system ( HHi) . In analogy with the procedure used in 
Section II, we find for each of the subalgebras listed in the previous section the associated 
five invariants along with the change of variable that must be substituted into system (HH!) 
in order to obtain the set of reduced equations. In each case, the symmetry variable is 
labelled by a. The invariants and change of variable are listed in Tables V and VI, while 
the system of reduced ordinary differential equations are listed in Tables VII and VIII. 
Where it is possible, we also give an explicit solution of the reduced system. We now 
proceed to describe the various classes and types of solutions which we obtain through 
solving the systems of reduced equations. 

We begin by listing those subalgebras for which we do not obtain an explicit solution 
in closed form. In the case of subalgebra £4, the fermionic fields ^ and must be constant 
and the bosonic field G must be related to F through the expression 

where F satisfies the differential equation 

(2(t^/3f - 3(tF2) F,, + (9(tF - Sa^'^) {F^f + {a^'^'F - 2F^) F, = 0. (54) 

Equation (15^ does not have the Painleve property. Similarity, we do not obtain an explicit 
solution in the case of subalgebra C21, which is closely related to the reduction obtained 
for subalgebra £4. 

For subalgebra Ci, we obtain the trivial solution where the four fields R, S, ^ and ip 
are constants. This is also one of the solutions of the reduced equations for subalgebra 

u 

A number of subalgebras lead to reductions which result in linear and travelling wave 
solutions in certain bosonic or fermionic fields. One of the solutions for subalgebra £5 
consists of a linear travelling wave in the fermionic field ^ of the form C,{x, t) = K(x — Cnt), 
while the other three fields, R, S and ip are constants. For subalgebra Cq, we obtain a 
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similar travelling wave in the field ip. Prom subalgebra £i8, we obtain the solution 
R{t) = r]iri2t + Di, S{t) = -r]ir]2t + D2, 



(55) 



^(x, t) = -qix + 7^1-02 ^ + -D3, ilj{x,t) ^ r}2X + r}2Dit + D4, 

where Di, Di, Di, Di are arbitrary fcrmionic constants. This solution is linear in t for 
the bosonic fields R and S and includes linear travelling waves in both fcrmionic fields. A 
similar solution is present for subalgebra £19. The solutions corresponding to subalgebras 
£22 and £23 are given by 

R = Ci, S = C2, ^{x,t) =r]2{x-C2t), iP{x,t)=r]2ln{x-Cit), (56) 

and 

R^Ci, S^C2, C{x,t)^rn^n{x-C2t), i^{x,t) ^ ru{x - C^t). (57) 

respectively. Here, the fcrmionic fields include both linear and nonlinear travelling waves, 
containing logarithmic functions. 

Another type of solution which we obtain quite frequently is the centered wave con- 
sisting of the ratio x/t. The solution corresponding to subalgebra £26 is given by 

i?(a;,t) = |, S{x,t) = j, i{x,t) = Q, 2P{x,t)=er]2^, (58) 

has simple centered waves in bosonic and fermionic fields. The solution corresponding to 
subalgebra £27 is very similar. Subalgebras £5 and Cq also provide us with the solutions 

R{x,t) — ^, S{x,t) — ^, ^ = 0, i^{x/t) an arbitrary function, (59) 



and 



x X 

R{x,t) = —, S{x,t) — —, ip — 0, Ci^/t) an arbitrary function, (60) 



respectively. Here, the bosonic fields R and S still consist of simple centered waves, but 
in each case the nonzero fermionic field consists of an arbitrary function of the ratio x/t. 

A number of solutions consist of functions similar to centered waves, but which 
contain important difi^crences. For example, subalgebra £j, in the case where a =7^ 
0, —1, —1/2, —1/3, leads us to the following solution 

4a^(l + 2a)^ (l+3a)(l+4a) (l+4a) 

4)(X t) — ^ LnX 2»(l+2a) t 2a + ^)^ 

(l + 3a)(l + a)2(12a3-2a2-5a-l)^ ^— ' 

15 



where Kq and Lq are defined in such a way that KqLq — 0. Here, the bosonic fields R 
and S still constitute centered waves, but the fermionic fields ^ and ip consist of specific 
powers of x and t which cannot be expressed in terms of the ratio x/t. This is also true of 
subalgebra Cg, for the case where a 0, —1, —1/2, —1/3, where we obtain the solution 

o/ fl + 2a\x fl + 2a\x 



1 + a J t' ^'^ \l + 3aj t' 

*' = (1 + 3a)(l + a)^(12a3 - 2a^ - 5a - i) j^^^^''^ + Rl. 

where /Cq and Lq are defined in such a way that LqKq = 0. For subalgebra £9, we obtain 
the solution 

-a(l + 3a) 2a a+3 2 

^{x,t) = Xqx 1-°'' ^1-°, i^jx^t) = Lpxi^t"-!. 

Other modified centered wave solutions include those for subalgebras £11, £12 and £13 
which are given by 

R{x,t) = - + e, S{x,t) = --e, ^{x,t) ^ Kt^^'^e^'^ , ^(x, t) = Li^/V^'^ 
R{x,t)^--e, S{x,t)^- + e, ^{x,t) = Kt^/^e-^J , iP{x,t) = Lt^/^e^'^ , 
and 



(64) 
(65) 



i?(x,t) = |-i£, S{x,t)^j + \s, ^{x,t) = Kt'/^e-'^, ^{x,t)^Lt'/^e'^, 

(66) 

respectively. 

Solving the reduced equations corresponding to subalgebra £3, we obtain the solution 

X C2 X C2 

R{x, t) — J + Cit + — , S{x, t) — — — Cit + — , ^ and ip are constants, (67) 

in which the bosonic fields R and S represent time-modified centered waves. The related 
subalgebra C20 gives us 

i?(a;, t) = J + Cit + , 5(a;,t) = J - Cit + , 
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in which time- modified centered waves are present in all fields (both bosonic and fermionic). 



For subalgebra £7, the case where a = —1/2 is handled separately. Here, the invariants 

of the subalgebra generator are 

X, tR, tS, t^^, -0, (69) 
and the required change of variable is 

R^\f{x), S^^-G{x), e = ^A(^), i^^i^ix). (70) 
This leads to the reduced equations 

-F + GF^ + V'^A^ = 0, -G + FG^ + A^V'^ = 0, -2A + GA^ = 0, = 0, 



(71) 



and we obtain the exponential solution 



R{x,t) = ^ S{x,t)- ^ ^ 



Ci ' ^ ' ^ C, ' (72) 

t) = Ko{l- e^^(^+^^))' , ^|Jix, t) = K,. 
A similar analysis for subalgebra £7, in the case where a = —1/3, leads to the exponential 
solution 

_ Goxe-^°(*-*") Cqx 

^{^^ ^) - l_ e-Co{t-to) ' ^K^^'^) - ^_ g-Co(i-«o) ' 

/ e-Co(t-to) \3 ^"^^^ 

In analogy with subalgebra £7, wc consider the equivalent cases for subalgebra Cg. The 
cases where a = —1/2 and a = —1/3 lead to the exponential solutions 

R{x,t)^^ 7^ ^, S{x,t)^^ ^, 

<--i Gi (74) 

ax, t) = ^, i;{x, t)^K^{l- e^^(-+^^))' , 



and 



ax,t)^Ko, ip{x,t)^Kix^ 



(75) 



1 _ e-Co(t-to) ' 1 _ g-C(,(t-to)' 

g-Co(t-t„) N 3 

respectively Exponential solutions are also obtained in the cases of subalgebras £24 and 
C25, for the case where a — —1/2 and are given by 

R{x, t) = ^ (e^^(^+^^) - 1) , S{x, t) = (e-^i(-H-^^) - 1) , 

Olt Oil 

1 n /I _ p-Cl(x + C2)\ C^S) 

ax,t) = ^!Z2 (1 - e^^(^+^^))' , V'(^,i) = 2^ In 
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and 

R{x, t) = --i- (e-^^(^+^^) - 1) , S{x, t) = -i- (e^^(-+^^) - l) , 

' ' (77\ 

respectively. 

Certain solutions are also obtained in terms of radicals. For subalgebra Cu, we obtain 
the solution 

2e 



R{x, t) = JCi (x- - let^) +C2 + et, S{x, t) = -yr\ Ci {x - let^) +C2 + et, 

V Ci V ^^g^ 

^ and are constant, 
and in the case of subalgebra £28, for a — —1/2, the following solution is present 

= (1 - + et, 

S(.. t) ^ -y^e. - - C„ (1 + ^p^;_3f _ + et, ^^^^ 

i{x, t) = Cir]iy/2ex - - Co + met + Di, 

ip{x,t) = --^rj2y/2sx - - Cp + r)2et + D2. 
Ci 

For subalgebra £10, we obtain an interesting solution to the reduced equations in the 
monomial form 

F{a) = /loa^ G{a) = Bqct', A(a) = Koa"\ n{a) = Loa^, (80) 
where the constants are determined to be 



LoKo 



l + a + 36' l + a + 36' l + a + 36' 

[p(l + a + 36) - 1] (1 + a + 26) 
p(l + a + 36)2 ' 
1 + a + 26\ / (1 + a + 36)p - (3 + 2a + 66) 



l + a + 36y \{l + a + 3b)p-{3 + 3a + 6b)^ 
(1 + a + 2bf{pa - 1 + 3p6 +p)(a2 - apfe + a + bab - pb - 3pb'^ + 66^ + 36) 



(1 + a + 36) V(-3 - 3a - 66 + p + pa + 3^6)^ 



(81) 



and p is the solution to the quadratic equation 

(-2a26 - 4a6 - 12a62 - 126^ - 186^ - 26)^^ 

+ {a^ + a^ + 10a26 + 33a62 + 12a6 - a + 366^ + 26 - 1 + 276^)p (82) 
+ (a^ + 4a + 5a6 + 66^ + 96 + 3) = 
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Finally, we note that the reduced equations corresponding to subalgebra £2, possess 
certain solutions expressed in terms of arbitrary functions. Specifically, we have the 
following three possibilities 

R = 0, ^ is a constant, while S{x) and 'ip{x) are arbitrary. (83) 

S = 0, i/) is a constant, while R{x) and ^(x) are arbitrary. (84) 

R = 0, S = 0, ^(x) and ipi^x) are bound by the condition C.x^'x = 0. (85) 

These functional degrees of freedom allow us to consider all kinds of interesting physical 
phenomena, including bumps, kinks and various types of multi-solitonic waves. 



VI Conclusion 



In this paper, we have successfully constructed a supersymmetric extension of the hydro- 
dynamical system in Riemann invariants (fTOj) . For both the standard and supersymmetric 
versions of (|TOl) . the Lie algebra of infinitesimal symmetries was computed and a system- 
atic classification of the subalgebra structure was performed. Through the use of the 
symmetry reduction method applied to both bosonic and fermionic variables, a number 
of interesting group-invariant solutions of both the classical and supersymmetric hydrody- 
namical models were determined. In the case of the supersymmetric model, the solutions 
included travelling waves, pure and modified centered waves, monomial, exponential and 
radical solutions and well as bumps, kinks and multi-solitonic solutions. The classical 
hydrodynamical model possesses many of the same types of solutions (time-modified cen- 
tered wave, exponential and radical solutions), but also solutions expressed implicitly. 
Conversely, it does not admit the solutions containing arbitrary functions which were 
present in the supersymmetric case. 

It is interesting to note that while the Lie symmetry algebra L of the classical hydro- 
dynamic model includes a generator J identified as a reverse boost, this symmetry is not 
present in the Lie superalgebra C of the supersymmetric model. Consequently, £ is a 
solvable algebra whereas L is not. 

The version of the symmetry reduction method presented in this paper has proved 
to be a useful tool for the purpose of constructing new, explicit, interesting solutions 
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of both the classical and supersymmetric versions of the hydrodyiiamical system in Rie- 
mann invariants. Further application of the methods used in this paper could lead to 
many physically interesting results for hydrodynamical systems, including the blow-up 
phenomenon. Particular analytic solutions have several advantages. For example, they 
allow us to observe qualitative behavior which might otherwise be difficult to find if only 
numerical descriptions are available. 
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Table I: Supercommutation table for the Lie superalgebra L spanned by the vector fields ([T8]) 





IVli 


iVi2 


w 


T 
J 




J-0 


Ml 














-Ti 


-To 


Ma 








-2W 


2J 


-Ti 


To 


W 





2W 





-Ma 





-Ti 


J 





-2J 


M2 





-To 





Ti 


Ti 


Ti 





To 








To 


To 


-To 
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Table II: Invariants of the one-dimensional subalgebras of L 



Subalgebra 


Invariants 


Relations and Change of Variable 


Li = {Ti} 


a = t, R, S 


E = R(t), S = S{t) 


L2 = {Ml} 


a=f,R,S 


R = R{a), S = S{a) 


L3 = {M2} 


a = xt, t^R, t^S 


R = ^F{a), S = ^G{a) 


L4 = {W} 


(J = t, R — J , S — J 


R = F{t) + 'j, S = G{t) + ^ 


L5 = {W- J} 


a = x'^ + 1^, 


T-) x tan 7^(f7) — t 
^ ~ x+tt&nF{cj) 1 




arctan R + arctan (|:) , 


Q xtaiiG(cF)—t 

'-' ~ x+tta,nG{a) 




arctan S + arctan , 




Le = {M2 + kMi,k^ 0,1,-1} 


a = xr^'i^', t~^T^>R, t~^k=iJS 


R = t^k^)F{a), S = t^T^)G{a) 


Lr = {Ml + M2} 


^ = t' 


R = xF{t), S = xG{t) 


Ls = {M2 - Ml} 


a = X, tR, tS 


R= \F{x), S = lG{x) 


Lg = {W + eMi,e = ±1} 


a = J — £ In t, R — elnt, S — elnt 


R = F{a) +elnt, S = G{a) + eliit 


Lio = {W - J + kMi,k ^ 0} 




7p XZ&Ylr{(J) — t 
~ x+tta.nF{cr) ' 




arctan R + arctan (|:) , 


^ X tan G{a)—t 

x+t tan G{(j) 




arctan S + arctan (|:) , 




Ln = {\V + sT,} 






L12 = {Ml + M2 + eTo} 


a = xe-^^\ \R, is 


R = xF{a), S = xG{a) 


Lis = {M2 - Ml + eTi} 


a = ie^^^, tR, tS 


R=}Fia),S = }Gia) 



24 



Table III: Reduced Equations obtained from the one-dimensional subalgebras of L 



Subalgebra 


Reduced Equation(s) 


Li = {Ti} 


i^t = 0, St = 


L2 = {Ml} 


-aRa + SRa = 0, -aSa + RS^ = 


L3 = {M2} 


-2F + aF^ + GF^ = 0, -2G + aG^ + FG^ = 


L4 = {W} 


Ft + ^G = 0, Gt + ^F = 


L5 = {W- J} 


-1 - (tan F)2 + 2cr(tan G)F^ + 2cr(tan G) (tan F)'^F^ = 0, 
-1 - {tanGf + 2a{tan F)Ga + 2a{tan F){tanGfGa = 


LQ = {M2 + kMi,k^ 0,1,-1} 


- |±iaF. + GF. = 0, ^G - |±iaG. + FG. = 


= {Ml + M2} 


Ft + GF = 0, Gt + FG = 


Lg = {M2 - Ml} 


-F + GF, = 0, -G + FG^ = 


Lg = {W + eM-i} 


-aF^ - eF^ + GF^ + £ = 0, -aG^ - eG^ + FG„ + e = 


Lw = {W -J + kMi,k^O} 


ka{secFfF^ + a taxiG {sec FfF^ - (tanF)^ -1 = 0, 
ka{secGfGa + crtanF(secG)2G^ - (tanG)^ -1 = 


Lu = {W + eTi} 


Ft + j^,G = 0, Gt + ^F = 


L12 = {Ml + M2 + eTo} 


-2eo-F^ + GF + aGF^ = 0, -2eaG„ + FG + ciFG,^ = 


Li3 = {M2 - Ml + eTi} 


-F + aF^ + 2£C7GF^ = 0, -G + aG^ + 2eaFGa = 
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Table IV: Supercommutation table for the Lie superalgebra H spanned by the vector fields ([19|) 





Ui 


U2 


U3 


r> 


Po 


Pi 


"V" 

11 


"V" 

12 


Di 











B 


-3Po 


-2Pi 








D2 














-Po 


-Pi 


-Yi 





D3 














-Po 


-Pi 





-Y2 


B 


-B 











-Pi 











Po 


3Po 


^0 


^0 


i^i 














Pi 


2Pi 


Pi 


Pi 

















Yi 





Yi 




















Y2 








Y2 
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Table V: Invariants of the one-dimensional splitting subalgebras of C 



Subalgebra 


Invariants 


Relations and Change of Variable 


A = {Pi} 


a = t, R, S, ^, 


R = R{t),s = s{t),c = m,i^ = m 


A = {Po} 


a = X, R, S, ^, ■0 


R = R{x), S = 5(x), ^ = ({x), ^ = il){x) 


A = {B} 


a = t,R-f,S-f,^,i^ 


R = Fit) + f,S = Git) + l, 
^ = i{t),^P = ^p{t) 


A = {Dx} 


a = ^,t'/^R, t^^S, 


R = t-^/^F{a), S = t-^/'^G{a), 
( = ^{a),i; = ^p{a) 


A = {^2} 




R = R{a), S = S{a), C = tA{a), ^p = ^l){a) 


A = {D^} 


a = J, R, S, ^, 1 


R = R{a), S = S{a), ^ = ({a), V' = tn{t) 


A = {-D2 + a£»i} 




R = tTff^F(cr), S = ti^G{a), 
^ = ti+3^A((7), tp = tp{a) 


A = + aDx} 




R = t^F{a), S = t^G(cr), 
^ = ^{a),i^ = t^n{a) 


A = {^3 + aD^} 


a = ^,R, S, t^^, 


R = R{a), S = S{a), 
^ = t^A{a), ip = ti+^ri(cr) 


>Cio = {Dz + aD2 
+6^1} 


, l+a+3fc fc 
6 -a -1 


-b -b 
R = ti+<^+3b i<\a), = ti+"+ai'G(o-), 

a 1 

^ = t^+'^+^b A[a), ip = t^+<^+'i'>Q,{a) 


Cxi = {^2 + eB} 


a = J — elnt, R — elnt, 
5 — elnt, |, 7/; 


R = F{a)+elnt, S = G{ct) + elnt, 
e = tA{a), V = i'ia) 


C12 = {D^ + eB} 


(T = | — elnt,i? — elnt, 
.5 -elnt, e, 1 


R = F{a)+elnt, S = G{a)+elnt, 
C = V = tn{a) 


jO.13 = {D3 + aD2 
+eB} 


0- = (1 + a)f - elnt, R - j^lnt, 
S-^lnt, t^C, t^V 


R = F{a) + ^ Int, S = G{a) + j^lnt, 
^ = t^A{a), tp = t^9,{a) 


jCu = {B + ePo} 


a = X — ^et^, R — et, S — et, ^, tp 


R = F{a) +et, S = G{a) + et. 


As = {1^1} 


X, t, R, S, tp 


N/A 


Ae = {Y2} 




N/A 


j0i7 = {Yi + eY2} 


X, t, R, S, ip — e^ 


N/A 
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Table VI: Invariants of the one-dimensional non-splitting subalgebras of L 



Subalgebra 


Invariants 


Relations and Change of Variable 




a = t, R, S, 
$, - rjix, i) - r]2X 


R = R{t), S = S{t), ^ = A{t) + riix, i/j = n{t) + r}2X 


>Ci9 = {Pq + ri^x 


a = X, R, S, 

C - vit, ^ - mt 


R = R{x), S = S{x), ^ = A{x) + r]it, ij; = n{x) + r]2t 


>C2o = {S+!?in 
+r/2>'2} 


a = t, R - J, S - J, 


R = F{t) + ^,S = Git) + ^, 
^ = A{t) + vi^,^ = Q{t) + mj 


+??2>2} 


a = ^,t'/^R, t'/^S, 
^ - Int, ij - ir/2 Int 


R = t-^/^Ria), S = t-^/^S{a), 
^ = A(cj) + ir/i Int, = n{a) + ir^alnt 


>C22 = {£'2 + ?y2>2} 


a = ^, R, S, |, V-??2lnt 


R = R{a), S = S{a), C = tA{a), = + r?2 Int 


i-23 = i-t^3 + 'yiJ^ij 


a = J, K, b, ^-rji Int, f 


K = K(u), b = byu), c, = A(^(7j -|- 771 Int, ip = t\l[a) 


-C24 = {-D2 + aDx 
+V2Y2} 


L + 'Aa 

a = t^R, t^S, 
tTTs-a^,^-r,2^^ 


R = t^F(o-), S = t^G{a), 
^ = tT+3^Aia),^ = n{a)+r!2^ 


jO-25 = {D3 + aDi 
+r?i>i} 


a= ,t^R, t^S, 


R = tT^F{a), S = t^G(c7), 
^ = A{a)+rn^^,i; = t^^n{a) 


C-2& = {D2 + eB 

+mY2} 


a = J — elnt, R — elnt, 

5- clnl, f, V - //2I11/ 


R = F{a)+elnt, S = G{a)+elnt, 

( = iMa), V = 9.{a) + //2I11/ 


C27 = {D3 + eB 
+V1Y1} 


a = J — elnt, R — elnt, 
S — elnt, ^ — rjilnt, j 


R = F{a)+elnt, S = G{a) + elnt, 
^ = A(cr) -1- 771 Int, V' = tn{a) 


£28 = {B + ePo 
+riiYi + ^212} 


a = X — \et^, R — et, S — et, 
^ - rjiet, - T]2et 


R = F{a) +et, S = G{a) + et, 
^ = A(cr) + riiet, tp = n{a) + ^et 
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Table VII: Reduced Equations obtained from the one-dimensional splitting subalgebras of L 



Subalgobra 


R (xiiiccxi Equation (s) 


A = {^'i} 


i?t = 0, St = Q, 6 = 0, Vt = 


A = {Po} 


5i?^+V'xCx = 0, + CxV-x = 0, SCx = 0, i?Vx = 


A = {B} 

•J L J 


Ft + iG = 0, Gt + iF = 0, Ct = 0, Vt = 


U = {Di} 


-\G - 2aGa + 3fT2/3FG(, + 9a^/3^<,Va = 0, 
-20-^^ + 3cj2/3G^^ = 0, -2ailj„ + Scr^/^F^^ = 


A = {D2] 


-aR„ + S'i?^ + VaA^ = 0, -crS*^ + RS„ + A^Va = 

A - CrA^ + S'A^ = 0, -Tlpa + -R-fAa = 


A = {D^} 


-aR^ + SR„ + Vl^i^ = 0, -(75^ + RS„ + ^a^^a = 
-aia + S'^^ = 0, n - (7fi(^ + = 


C7 = {D2 + aDi} 


-j^F - (1 + 2a)aF^ + (1 + ^a)a'^^)GF^ + (1 + 3a)V(T7^)VaA^ = 0, 
-j^G - (1 + 2a)c7G^ + (1 + 3a)c7(^)FG^ + (1 + 3a) V(^)a^V;^ = 0, 

1 3a 

- (1 + 2a)(7A<, + (1 + ^a)a^Gka = 0, 

3a 

-(1 + 2a)aip^ + (1 + 3a)(7T+3^FV'a = 0, 


jCs = {D3 + aDi} 


-j^F - (1 + 2a)(7F^ + (1 + 3a)a(TT^)GF<, + (1 + ^af a^^) VL^^^ = 0, 
- ii^G - (1 + 2a)aGa + (1 + 3a)(T(^)FG<, + (1 + ?,af a^W) ^^Q.^ = 0, 
-(1 + 2a)ai^ + (1 + 3a)f7i+5^Ge<x = 0, 

T 3a 

yq^3^fi - (1 + 2a)crO^ + (1 + ^a)a^F^a = 0, 


Lg = {D3 + aD2} 


-Cri?cr + SRa + ^aK = 0, -(tS^ + i?^,^ + K^a = 

^A - aK + SK = 0, - (J^a + Rn^ = 


jCio = {D3 + aD2 
+bDi} 


-TT^^ - ^F'^ + {mi) + {mi) ^^(™^)iuv. = 0, 


-T+^G - aG, + (M) ^(i+.+3Ji.'G. + {^) a(i+»+-)A.n. = 0, 

/ \ / b \ 

TT^A - + (mi) ^^^+--^GA. = 0, 


1+0+36" ' 1^1+0+26^" '-TliCT — U, 


-Cii = W2 + eB} 


-(tF^ - eFa- + GFa + e + tha^a = 0, 
-crG^ - £G(, + FGa +e + K'^a = 0, 
A - (tA^ - eAa- + GAa^ = 0, -(ripa - ^i^a + Flpa = 


jCi2 = {D3 + sB} 


-aF^ - eF„ + GF^ + e + ^^^a = 0, 

-aG^ - sGa + FG^ + e + ^Ma = 0, 
-(T^a - eia + = 0, 17 - ctO^ - eJl,, + Ff]^ = 


jC-13 = {D3 + aD2 
+eB} 


-aF^ - eF, + (1 + a)GF, + ^ + (1 + a)2j],A, = 0, 
-uG^ - sGa + (1 + a)FG„ + + + afKVL^ = 0, 
^A - aK - eAa + (1 + a)GA^ = 0, 
j^n - an„ - efl^ + (1 + a)Fn^ = 


Cii = {B + ePo] 


GF„ + e + tpa^a = 0, FG^ + e + i^^^ = 0, 
G^^ = 0, FV-a = 
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Table VIII: Reduced Equations obtained from the one-dimensional non-splitting subalgebras of 



Subalgcbra 



Reduced Equation(s) 



>Ci8 = {Pi + mYi + mY2} 



Rt + mm 

At + Srii 



0, 
0, 



St + Vim 
nt + Rm 



0, 




jO-19 = {Po + mYi + V2Y2} 



SRx + ^x-Ax — 0, 

m + SAx = 0, 



RSx + Ax^x — 0, 
r]2 + Rnx = 



C20 = {B + rjiYi + mY2} 



t^Ft + tG + mm = 0. 

tAt + mG = 0, 



fGt + tF 



mm 



mt + rj2F = 



C21 = {Di + mYi + mY2] 



-\F- laF„ + GF„ + n^A^ = 0, 
- laG^ + FG^ + A^O^ = 0, 
■laA„ + \m + GA^ = Q, -laVL^ + \m + FVL„ = Q 



C22 = {D2 + mY2} 



-aRa + SRa + n^Aa- = 0, 
A - crAo- -I- SA^ = 0, 



-aSa + RSa + A^n^ = 

-a^a + m+ R^a = 



j0.23 = {D3+mYl} 



-aRa + SRa + ^aAa = 0, 

-aAcr + m + SA^ = 0, 



-aS^ + RSa + A^n^ = 

Q - an^ + Rn^ = 



JO-24 = {D2 + aDi + 772 12} 



AF-^F'r + ]^,a(TTrJGF,+ [^) a(T^)o.A. = 0, 

T^aG - + i±ia(T^)FG. + l^yaiT^^)A,n, = 0, 
1 A - aA„ + ii^cr i^GA„ = 

l+3a " <^ ^ l+2a <^ ' 



-crn, + j^m + ^aTfs-aF^}, = 0, 



C25 = {D3 + aDi+mYi} 



l+3a 



k-aG - ^G. + l±f|a(i^)FG. + (0^0^^) A,^, = 0, 

0, 



~ctA^ + 



1 71, + l+3ao-l+3^GA 
l+3al/l ^ l+2a" ""^^^f^ 

l+3a , 



l+3a 



» - + = 0, 



C26 = {D2 +eB + mY2] 



-aFa - sF^ + GFa + s + n^A„ = 0, 
-aGa - sGa + FGa + £ + A^^ = 0, 

A - rxA, - cA, + GA, = 0, -a9., - cQ, + FO_^ + //2 = 



>C27 = {D3 + eB + mYi} 



-aAn 



-aF^ - eFa + GF^ + e + = 0, 
-aGa - eGa + FG^ + e + A^^l^ = 0, 
eAo- -I- GA^ + r/i = 0, Vt- a^l^ - eQ„ + F^^ = 



£28 = {B + ePo + r/iYi + 7/2^2} 



GP<, + £ + O^A,, = 0, 
GA„ + m£ = 0, 



FG^ + e + A^n^ = 

Fn^ + m£ = o 
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